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Abstract 

Peristaltic transport of copper-water nanofluid in an inclined channel is reported in the presence of mixed convection. Both 
velocity and thermal slip conditions are considered. Mathematical modelling has been carried out using the long 
wavelength and low Reynolds number approximations. Resulting coupled system of equations is solved numerically. 
Quantities of interest are analyzed through graphs. Numerical values of heat transfer rate at the wall for different parameters 
are obtained and examined. Results showed that addition of copper nanoparticles reduces the pressure gradient, axial 
velocity at the center of channel, trapping and temperature. Velocity slip parameter has a decreasing effect on the velocity 
near the center of channel. Temperature of nanofluid increases with increase in the Grashoff number and channel 
inclination angle. It is further concluded that the heat transfer rate at the wall increases considerably in the presence of 
copper nanoparticles. 



Citation: Abbasi FM, Hayat T, Ahmad B, Chen G-Q (2014) Slip Effects on Mixed Convective Peristaltic Transport of Copper-Water Nanofluid in an Inclined 
Channel. PLoS ONE 9(8): e1 05440. doi:1 0.1 371/journal.pone.01 05440 

Editor: Victor M. Ugaz, Texas A&M University, United States of America 

Received February 15, 2014; Accepted June 30, 2014; Published August 29, 2014 

Copyright: © 2014 Abbasi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits 
unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 

Funding: This paper was funded by the Deanship of Scientific Research (DSR), King Abdulaziz University, Jeddah, under grant no. 9-130/1434 HiCi. The authors, 
therefore, acknowledge with thanks DSR technical and financial support. The funders had no role in study design, data collection and analysis, decision to publish, 
or preparation of the manuscript. 

Competing Interests: The authors have declared that no competing interests exist. 
* Email: abbasisarkar@gmail.com 



Introduction 

Mixed convective flows have received numerous attention from 
the researchers the world over due to their wide industrial and 
engineering applications. Convection is seen in the ocean currents, 
sea-wind formation, rising of plume of hot air from fire, formation 
of micro-structures during the cooling of molten metals, solar 
ponds and in fluid flows around heat dissipation fins. Convection 
in channels plays a vital role in heat exchangers, removal of 
nuclear waste and in modern cooling/heating systems. "Nano- 
fluids" refer to heat transfer liquids with enhanced heat transfer 
capability. Such materials are obtained by suspending nanopar- 
ticles in traditional heat transfer liquids e.g. water, ethylene glycol 
and engine oil. A number of experimental studies elaborate that 
the convective heat transfer in nanofluids is considerably higher 
than that of the base fluids. Enhanced heat transfer in nanofluids 
has enabled there use in several electrical and engineering 
applications. This gives rise to a new branch of mechanics named 
nanofluid mechanics. The wide utility of nanofluids is the reason 
for growing interest in nanofluid mechanics by the researchers of 
modern era. The term nanofluid was introduced by Choi [1] in his 
effort to enhance the thermal conductivity of the fluids. 
Buongiorno [2] proposed that thermophoresis and Brownian 
motion are significant in the dynamics of nanofluids. Some studies 
using the Buongiorno's model of nanofluids can be seen through 
refs. [3-5]. Expressions for the effective density and viscosity of the 
nanofluids were computed by Xuan and Li [6]. Brinkman's [7] 
model for the effective viscosity of the two phase flow was found to 
be in good agreement with the experimental results (Xuan and Li. 
[8]). Tiwari and Das [9] proposed a model for the analysis of 



nanofluids employing the results of Xuan and Li [6] and Brinkman 
[7] . Model of Tiwari and Das [9] was used by many researchers 
for the analysis of nanofluid dynamics and was found in good 
agreement with the experimental data (see [10-13]). 

Peristalsis is a well-known mechanism for fluid transport in 
physiology. In this mechanism sinusoidal waves travel on the walls 
of tube like organs of human beings propelling the fluid contained 
within the tube in the direction of their propagation. In physiology 
the principle of peristalsis is seen in the transport of food through 
oesophagus, movement of chyme in intestines, urine transport 
from kidneys to the bladder, bile transport in bile duct, transport of 
spermatozoa, vasomotion of blood vessels and many others. 
Peristalsis has proven very useful in fluid transport over short 
distance preventing the fluid from being contaminated. Subject to 
its utility, such mechanism has been adopted by the engineers in 
designing several industrial appliances including roller and finger 
pumps and peristaltic pumps in heart lung and dialysis machines. 
Latest hose pumps of several kinds operate through the principle of 
peristalsis. Peristaltic type flow is readily being used in the nuclear 
industry for the transport of corrosive fluids. Subject to such wide 
occurrence of peristalsis in physiology and industry many 
theoretical investigations are carried out in different flow 
configuration [14-22]. It is now well admitted fact that all the 
tubular organs facilitating fluid flow in human body are internally 
lubricated with mucus and secretion layers. These layers in turn 
prevent the fluid from sticking to the walls. In such cases the no- 
slip condition between the fluid and solid boundary is not valid. 
With such motivation some authors investigated the peristaltic 
transport of traditional fluids with slip effects [23—25]. 
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Figure 1. Geometry of the problem. 

doi:1 0.1 371 /journal. pone.01 05440.g001 

It is noticed from existing information on the topic that not 
much has been said yet about the peristaltic motion of nanofluids. 
Few recent investigations constructed such flow models (see [26— 
27]). In these attempts, the Buongiorno's model of nanomaterial 
has been employed. However there is not a single attempt which 



can address the peristaltic transport of nanofluid employing Tiwari 
and Das [9] model. Intention in present communication is to 
prepare such fluid model describing peristalsis. Hence we consider 
here the copper-water nanofluid to discuss the peristaltic flow in a 
symmetric channel with slip conditions. The channel walls satisfy 
the velocity and thermal slip conditions. Inclined channel filled 
with an incompressible fluid is considered. Relevant modeling is 
presented. Computations for the governing problems are made 
numerically. Plots for different flow quantities of interest are 
displayed and analyzed. 

Formulation and Solution 

Consider the peristaltic flow of copper- water nanofluid in a 
symmetric channel of width 2d. The channel is inclined at an 
angle a to the horizontal. Moreover the channel walls are 
maintained at constant temperature To. Cartesian coordinates are 
chosen in such a way that the X-axis lies along the length of the 
channel whereas Y -axis is normal to X -axis. Waves of 
wavelength X and amplitude a\ travel on the channel walls with 
constant speed c. These waves are responsible for flow generation 
in the channel. Schematic diagram of the problem is given in 
Fig. 1. 

Wave geometry is of the form 



:H(X,1)- 



(2n — 
-_ci\ cos[ — (X- 



ct 



(1) 



where -\-H and —H denote the upper and lower walls 
respectively. Velocity profile for such a flow is 
\=[U(X,Y,t),V(X,Y,t),0]. The two-dimensional continuity 
equation for an incompressible flow can be written as follows: 



Table 1. Nomenclature. 




Nomenclature 


H (X,X) 


Peristaltic wall 


d 


Half channel width 


a n 


Amplitude of peristaltic wave 


X 


Wavelength of peristaltic wave 


U 


X-component of velocity 


V 


F-component of velocity 


P 


Pressure 


Pnf 


Density of nanofluid 


Veff 


Effective viscosity 


9 


Acceleration due to gravity 


Pnf 


Thermal expansion coefficient of nanofluid 


y. 


Channel inclination angle 


T 


Temperature 


To 


Temperature of the walls 


C nf 


Specific heat of nanofluid 


K e ff 


Effective thermal conductivity 


0> 


Dimensional heat source/sink parameter 


Hf 


Viscosity of base fluid 


Pf 


Density of base fluid 


Pp 


Density of nanoparticles 


c f 


Specific heat of base fluid 


c P 


Specific heat of nanoparticles 


Pf 


Thermal expansion coefficient of base fluid 


9 


Nanoparticles volume fraction 


Pp 


Thermal expansion coefficient of nanoparticle 


K f 


Thermal conductivity of base fluid 




Thermal conductivity of nanoparticles 


E 


Eckert number 


Re 


Reynolds number 


Pr 


Prandtl number 


Gr 


Grashoff number 


Br 


Brinkman number 


•A 


Stream function 


P 


Dimensionless velocity slip parameter 


F 


Dimensionless volume flow rate in wave frame 


y 


Dimensionless thermal slip parameter 


n 


Dimensionless volume flow rate in fixed frame 




Pressure rise per wavelength 


doi:1 0.1 371/journal.pone.01 05440.t001 
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Table 2. Numerical values of the thermophysical properties. 





Property 


Basefluid/water 


Nanoparticles/copper (Cu) 


Density (kg/m 3 ) 


997.1 


8933 


Thermal conductivity (W/mK) 


0.613 


400 


Specific heat (J/kgK) 


4179 


385 


Thermal expansion coefficient (1/K)x10 6 


210 


16.65 



doi:1 0.1 371/journal.pone.01 05440.t002 
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Figure 2. Pressure gradient for variation in nanoparticle volume fraction when Gr = 5, a = 0.7, 17 = 1.2, Br = 0.3, £ = 1.0, /? = 0.1, y = 0.1 and 
a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g002 
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Figure 3. Pressure gradient for variation in channel inclination angle when Gr = 5, « = 0.7, tj = 1.2, l?r = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 and 
0 = 0.1. 

doi:1 0.1 371 /journal.pone.01 05440.g003 
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Figure 4. Pressure gradient for variation in Grashoff number when 0 = 0.1, « = 0.7, tj = 1.2, Br = 0.3, £ = 1.0, /? = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g004 



dU dV fd e T7 d\ T7 dP (d 2 V d 2 v\ 

^ ^ ' () ^H^ + ^ + F ^J F= -^ + ^^ + ^J (4 ) 

where U and V are the longitudinal and transverse components —g^pfy ^(T—Tq) cos a. 
of velocity in the laboratory frame. The X and Y components of 

momentum equation in presence of mixed convection are Energy equation with viscous dissipation and heat generation/ 



(d — d - d \- dP d 2 U d 2 U\ 



absorption is given by 




Figure 5. Pressure gradient for variation in velocity slip parameter when Gr = 5, « = 0.7, j/ = 1.2, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 and 
a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g005 
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Figure 6. Pressure rise per wavelength for change in nanoparticle volume fraction when Gr = 5, a = 0.7, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 
and a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g006 




In these equations P is the pressure in laboratory frame, fi e ffthe 
effective viscosity of nanofluid, p n f the density of nanofluid, fi n f 
thermal expansion coefficient of the nanofluid, g the acceleration 
due to gravity, T the fluid temperature, C n f the specific heat of 
nanofluid, K e /f the effective thermal conductivity of nanofluid and 
0 the dimensional heat source/ sink parameter. Nomenclature of 
all the quantities is also provided through Table 1. Effective 
viscosity of the nanofluid can be computed with the help of 




Figure 7. Pressure rise per wavelength for change in channel inclination when Gr = 5, « = 0.7, Br = 03, £=1.0, 0=0.1, y = 0.1 and 0 = 0.1. 
doi:1 0.1 371 /journal.pone.01 05440.g007 
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Figure 8. Pressure rise per wavelength for change in Grashoff number when 0 = 0.1, a = 0.7, Br = 0.3, £ = 1.0, /? = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g008 



equations for two-phase flow. In this study we have used the 
relation proposed by Brinkman [7] 



2.5 " 



Relations for effective density, specific heat and thermal 
conductivity of nanofluid by Tiwari and Das [9] are 



p nf = (1 - #) Pf + <l>p p , (pC) nf = ((1 - <t>){Cp) f + <t>{Cp) p \ 



(6) (j>P) nf = (1 - 4>)PfPf + WpP p , 



K eff = Kp + lKf-l^Kf-Kp) 0) 
K f Kp + lKf + tiKf-Kp) ' 



Here Maxwell-Gamett's (MG-model) is used to approximate the 
effective thermal conductivity of the nanofluid. In these equations 
(f) denotes the nanoparticle volume fraction, py the density of base 




Figure 9. Pressure rise per wavelength for change in velocity slip parameter when Gr = 5, « = 0.7, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 and 

C£ = 7r/4. 

doi:1 0.1 371 /journal.pone.01 05440.g009 
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Figure 10. Axial velocity for different nanoparticle volume fraction when Gr = 5, « = 0.7, x = l, j/ = 1.5, Br = 0.3, £ = 1.0, /? = 0.1, y = 0.1 and 
a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g01 0 



fluid (water in this case), p p the density of nanoparticles, /?y the 
thermal expansion coefficient of base fluid, fi p the thermal 
expansion of nanoparticles, fij- the viscosity of base fluid, K p the 
thermal conductivity of nanoparticles and Kf the thermal 
conductivity of base fluid. Numerical values of these properties 
are readily available which we include in the Table 2 [9,10]. 

The quantities from the laboratory frame (X, Y, 1) to wave 
frame (x, y) are related by the following transformations: 



x = X-ctJ= Y,u = U-c,v= V,p(x,y) = P(X, Y,t), (8) 

in which w, v and p are the velocity components and pressure in 
wave frame (x,y). Equations (2)— (5) through Eqs. (6-8) yield 




Figure 11. Axial velocity for different Grashoff numbers when 0 = 0.1, a = 0.7, x=l, // = 1.5, Br = 0.3, £ = 1.0, p = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g01 1 
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Figure 12. Axial velocity for different inclination angles when Gr = 5, « = 0.7, x = l, q=1.5, Br = 0.3, £ = 1.0, ft = 0.1, y = 0.1 and 0 = 0.1. 
doi:1 0.1 371 /journal.pone.01 05440.g01 2 
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Figure 13. Axial velocity for different values of velocity slip parameter when Gr = 5, a = 0.7, x = l, ij = 1.5, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 
and a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g01 3 
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Figure 14. Behavior of streamlines for variation in nanoparticle volume fraction when Gr = 5, a = 0.7, rj = 1.5, l?r = 0.3, £ = 1.0, /? = 0.1, 

y = 0.1 and cc = 7r/4. 

doi:10.1371/journal.pone.0105440.g014 



((l-«P/ + ^)[ (w+c)^+v-^j(w+c) 



5/7 



d 2 w d 2 w 



(l_^) 25 ^x 2 dj 2 
+ g{Q-#)pfPf + WpP P ) (T-Tq) sin a, 



(( 1 - <t>) Pf + + 0 A + p A) v 



3/> 



e 2 v e 2 v 



5j (i-<^) 2 - 5 ^x 2 ej 2 

?(Q-4)PfPf + <l>PpPp) (T-T 0 ) cos a, 



(10) 



(11) 



((l-^XCp) / + «Cp) J ,)l(«+c)^+v4]7' 



dx dy_ 



d 2 T d 2 T 



+ 



(W) 



2.5 



+ 



'5m dvV 



With the help of following dimensionless quantities 



(12) 




Figure 15. Behavior of streamlines for variation in channel inclination when Gr=5, « = 0.7, // = 1.5, 5r = 0.3, £ = 1.0, /? = 0.1, y = 0.1 and 
0 = 0.1. 

doi:1 0.1 371 /journal.pone.01 05440.g01 5 
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Figure 16. Behavior of streamlines for variation in velocity slip parameter when Gr = 5, a = 0.7, 17 = 1.5, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 
and a = n/4. 

doi:1 0.1 371 /journal.pone.01 05440.g01 6 



x y u v d H a\ 
X d ,U c ,V cd' V d ,<2 d 



proximations Eqs. (10-12) are reduced to the following expres- 
sions: 



ck\if 



d 2 p T-T 0 /i f p f cd 
6 = — — — , i) = — , Re = — — , 



Ct 17 c o 
, E = — , Pr 



r o Pf y-f 
2 ^ Gr ^pjW^f_ 



(13) 



X ' CfTo 



-f 



H f c 



— = ^T^r-r + Gr sm oc < \ —<b + <b\ - — - 

dx (l-^) 25 ^ 3 \ ^ UP/0/ 



(14) 



(15) 



d 2 0 dx// dx// 



-,u= — ,V— , 
and applying the long wavelength and low Reynolds number ap- 



K p + 2K f -2<t>(K f -K p ) d 2 d Br d 2 A 
K p +2K f + t(K f -K p ) dy> + (1-^)" \ty ) +& 1 ° j 

where the continuity equation is identically satisfied and Eq. (15) 



2.0 



1.5 



CD 



1.0 



0.5 



















0=0.04, 


0.08, 0.12, C 


" * " - . *'x\ 

. X X 






% s • X 

s X \\ 

v xv 

N * • 








V 
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0.5 



1.0 



1.5 



Figure 17. Temperature profile for different nanoparticle volume fractions when Gr=5, a = 0.7, x = \, J/ = 1.5, Br =0.3, £ = 1.0, p = 0.1, 

y = 0.1 and a = n/4. 

doi:1 0.1 371/journal.pone.01 05440.g01 7 



PLOS ONE | www.plosone.org 



10 



August 2014 | Volume 9 | Issue 8 | e105440 



Peristaltic Transport of Copper-Water Nanofluid 




Figure 18. Temperature profile for variation in Grashoff number when 0 = 0.1, a = 0.7, x = l, // = 1.5, Br = 0.3, £ = 1.0, p = 0.1, 7 = 0.1 and 

C£ = 7r/4. 

doi:1 0.1 371 /journal.pone.01 05440.g01 8 



shows that p^piy)- In these equations ij/ is the stream function, 
Re the Reynolds number, Br the Brinkman number, E the Eckret 
number, Pr the Prandtl number, 3 the wave number, 6 the 
dimensionless temperature and e the dimensionless heat source/ 
sink parameter. From Eqs. 14 and 15 one obtains 



0 ^(r^v +Grsm T^ + ^PvJ|^- (17) 

Defining 77 and F as the dimensionless mean flows in laboratory 
and wave frames by 




Figure 19. Temperature profile for variation in channel inclination angle when Gr = 5, a = 0.7, x = l, // = 1.5, Br = 0.3, £=1.0, /? = 0.1, 
7 = 0.1 and 0 = 0.1. 

doi:1 0.1 371 /journal.pone.01 05440.g01 9 
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Figure 20. Temperature profile for variation in velocity slip parameter when Gr = 5, « = 0.7, x = l, tj = 1.5, Br = 0.3, £ = 1.0, 0 = 0.1, y = 0.1 
and ce = 7r/4. 

doi:1 0.1 371 /journal.pone.01 05440.g020 




Figure 21. Temperature profile for variation in thermal slip parameter when Gr = 5, a = 0.7, x = l, // = 1.5, Br = 0.3, £ = 1.0, ft = 0.1, 0 = 0.1 
and ce = 7r/4. 

doi:1 0.1 371 /journal.pone.01 05440.g021 
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Figure 22. Change in heat transfer rate at the wall for variation in nanoparticle volume fraction Gr=5, a = 0.7, v= 1. //= 1.5. Br = 0.3, 
£ = 1.0, /? = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g022 



= Q F= I_ 

cd ' cd 

two flow rates can be related as follows: 



rj=F+l, 



where 



.0 



dy. 



(18) 



(19) 



(20) 



The velocity and thermal slip conditions are 



U- U w = n/3 l T xy , T-T w 



■ ny l 



dT 

dy 9 



where U is the velocity of fluid, x xy the component of extra stress 
tensor, T the temperature of fluid, /? l9 yj are respectively the 
dimensional velocity and thermal slip parameters and n the unit 
normal on the wall. The dimensionless velocity and thermal slip 
conditions finally become 




Figure 23. Comparison of heat transfer rate at the wall of water and Cu-water nanofluid for change in heat source/sink parameter 
when Gr = 5, « = 0.7, x = l, t/ = 1.5, Br = 03, /? = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g023 



PLOS ONE | www.plosone.org 



13 



August 2014 | Volume 9 | Issue 8 | e105440 



Peristaltic Transport of Copper-Water Nanofluid 




Figure 24. Comparison of heat transfer rate at the wall of water and Cu-water nanofluid for change in Grashoff number when 

a = 0.7, x = l, J7 = 1.5, Bv = 03, £ = 1.0, /? = 0.1, y = 0.1 and a = n/4. 
doi:1 0.1 371 /journal.pone.01 05440.g024 



d 2 x/j 



30 
Ty 



where /?( = — ) and = are respectively the dimensionless 

velocity and thermal slip pramaeters. The appropriate boundary 
conditions for the present flow model are 



(A = 0,iA =0,^ = 0, ^ j/ = 0, 



(i-# 5 ^ 2 



50 (21) 
-1,0 + y^r- =0,aty = h. 
dy 



with A= 1 +acos(27ix). Pressure rise per wavelength (A/?;J is 
defined as follows: 



o uX 




0.0 0.1 0.2 0.3 0.4 0.5 



Br 

Figure 25. Comparison of heat transfer rate at the wall of water and Cu-water nanofluid for change in Brinkman number when 

Gr = 5, fl = 0.7, x = l, J7 = 1.5, £ = 1.0, /? = 0.1, y = 0.1 and « = tt/4. 
doi:1 0.1 371 /journal.pone.01 05440.g025 
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Table 3. Numerical values of heat transfer rate at the wall for change in the values of different parameters when Gr = 5, a = 0.7, 
x=l, 17 = 1.5, Br = 0.3 and y = 0A. 







a 


£ 


^ flYM 
~V () 


0.0 


7l/4 


1.0 


2.93204 


0.05 






2.79028 


0.10 






2.67728 


0.15 






2.58574 


0.10 


7l/6 




2.46308 




71/4 




2.67727 




n/3 




2.76587 




Till 




2.86097 




7l/4 


0.0 


0.78554 






1.0 


2.67227 






2.0 


4.70860 






3.0 


4.90948 



doi:1 0.1 371/journal.pone.01 05440.t003 



The system of Eqs. (14-16) subject to boundary conditions (21) 
is solved numerically using NDSolve in Mathematica. We have 
taken the step size equal to 0.01 for the variations in both variables 
x and y. Obtained numerical results are analyzed graphically in 
the next section. 

Results and Discussion 

Pressure gradient, pressure rise per wavelength, axial velocity 
and trapping are important quantities in the study of peristalsis. 
On the other hand the analysis of heat transfer phenomena is 
significant in the study of nanofluids. Hence the present section is 
prepared to analyze the behavior of pressure gradient, pressure 
rise per wavelength, axial velocity, streamlines, temperature and 
heat transfer rate at the wall for variation in different embedded 
parameters. 

Pressure gradient and pressure rise pre wavelength 

Pressure gradient is studied through Figs. 2-5. These Figs, show 
that pressure gradient is minimum in the occluded region of 
channel and attains maximum value in the wider part. Addition of 
copper nanoparticles reduces the pressure gradient. In addition an 
increase in the nanoparticle volume fraction decreases the pressure 
gradient. Such decrease is large in both the wider and occluded 
regions (see Fig. 2). Fig. 3 depicts that there is an increase in 
pressure gradient near the wider part of channel when inclinations 
increases. Effect of Grashoff number on pressure gradient is 
similar to that of channel inclination angle (see Fig. 4). Pressure 
gradient enhances when the velocity slip parameter is increased 
(see Fig. 5). Such increase is large in the occluded portion of the 
channel. 



Pressure rise per wavelength versus flow rate is shown in the 
Figs. 6-9. One common observation from these Figs, is that 
pressure rise per wavelength decreases with an increase in the flow 
rate. Fig. 6 shows that addition of nanoparticles decreases the free 
pumping flux (value of rj for Ap; i =0). Moreover an increase in the 
value of nanoparticle volume fraction enhances the pressure rise in 
retrograde pumping region (rj<0, A/?; >0). However reverse 
situation is seen in the augumented pumping region (rj>0, 
Apx <0). Figs. 7 and 8 show that increase in the values of Grashoff 
number and channel inclination angle result in the increase of 
pressure rise per wavelength. Fig. 9 shows that large velocity slip 
parameter decreases the peristaltic pumping region (fj>0, 
Ap;>0). Value of pressure rise per wavelength decreases largely 
in the retrograde pumping region when /? increases. 

Axial velocity and trapping phenomena 

Effects of several flow parameters on axial velocity are shown in 
Figs. 10-13. These Figs, show that velocity profile traces a 
parabolic trajectory with maximum value at the center of the 
channel. Fig. 10 depicts that nanofluid with low concentration of 
nanoparticles possesses higher value of the velocity at the center. 
Which concluded that higher nanoparticles volume fraction 
provides more resistance to the flow. However near the channel 
walls, this observation is not true which primarily is due to slip 
effects. Figs. 1 1 and 12 show that increase in Grashoff number and 
channel inclination angle result in an increase in velocity at the 
center. As in previous case this observation is also reversed near 
the channel walls. Velocity near the channel walls increases with 
an increase in the velocity slip parameter (see Fig. 13). However 
maximum value of velocity at the center decreases subject to an 
increase in ft. 



Table 4. Comparison of present study with the results obtained by obtained by Ali et al [28] when a = 0.6 and /? = 0.02. 

Alit et al. [28] Present study for Gr = 0, <j> = 0 

Critical value of v\ f or tsp k 0.45022 0.45233405 

Critical value of r\ for F k 0.3501 0.3505850 

doi:1 0.1 371/journal.pone.01 05440.t004 
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Trapping for different involved parameters is studied through 
Figs. 14-16. Some volume of fluid during flow tends to get trapped 
within streamlines. This phenomenon is usually known as trapping 
and the volume of fluid thus trapped is named as bolus. Size of 
such bolus reduces with an increase in the value of nanoparticle 
volume fraction (see Fig. 14). Hence it can be stated that addition 
of nanoparticles reduces the trapping. Size of the trapped bolus 
slightly increases with an increase in the channel inclination angle 
(see Fig. 15). Fig. 16 depicts that size of the trapped bolus reduces 
with an increase in the velocity slip parameter. 

Heat transfer analysis 

Impact of several embedded parameters on temperature and 
heat transfer rate at the wall are examined through Figs. 17-25. It 
is clear from Fig. 1 7 that an increase in the nanoparticle volume 
fraction shows a decrease in fluid temperature. This clearly 
witnesses that the heat being generated by the heat source is 
quickly transferred to the walls (due to addition of more 
nanoparticles) and consequently the fluid temperature decreases. 
High thermal conductivity of nanoparticles plays a key role in 
quick dissipation of the fluid temperature. This justifies the use of 
copper nanoparticles in different types as the coolants. Figs. 18 
and 1 9 show that increase in the values of Grashoff number and 
channel inclination enhance the fluid temperature. This is due to 
the fact that increase in Gr and a facilitates the convective heat 
transfer within the fluid. Figs. 20 and 21 show that velocity and 
thermal slip parameters have opposite effects on temperature of 
cu-water nanofluid. Temperature of nanofluid decreases with an 
increase in the velocity slip parameter. 

Heat transfer rate at the wall for variation in nanoparticle 
volume fraction is analyzed in Fig. 22. It is observed that larger 
nanoparticle volume fraction increases the heat transfer rate at the 
wall. Hence addition of copper nanoparticles facilitates the heat 
transfer between the fluid and solid boundary and ultimately 
preventing the fluid from heating. Figs. 23-25 present a compar- 
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ison of the heat transfer rate at the wall of water in the absence and 
presence of copper nanoparticles. All these Figs, show that copper- 
water nanofluid possesses higher heat transfer rate at the wall 
when compared with ordinary water. Also heat transfer rate is 
seen to increase by increasing 8, Gr and Br. Numerical values of 
heat transfer rate at the wall for variations in a and £ are given 
in Table 3. It is observed that heat transfer rate at the wall 
increases when a and e are increased. However it decreases 
through larger Comparison between the limiting case of present 
study and the analysis performed by Ali et al. [28] is provided 
through Table 4. Critical values of flow rate are provided and 
compared. Good comparison is found between the two. 

Conclusions 

Velocity and thermal slip effects on the peristaltic transport of 
copper-water nanofluid in an inclined channel are studied. 
Findings of present analysis indicate that the addition of copper 
nanoparticles reduces the pressure gradient. Increase in the value 
of nanoparticle volume fraction enhances the pressure rise in 
retrograde pumping region. Nanofluid with low concentration of 
copper nanoparticles possesses higher value of velocity at the 
center of channel. Velocity slip parameter has a decreasing effect 
on the maximum value of axial velocity. Trapping is reduced in 
presence of nanoparticles. Grashoff number and channel inclina- 
tion angle have increasing effect on the temperature of nanofluid. 
Increase in nanoparticle volume fraction largely enhances the heat 
transfer rate at the wall. 
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